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The baryon-baryon continuum invariant mass spectrum generated from 

relativistic nucleus + nucleus collision data may reveal the existence of 

doubly-strange dibaryons not stable against strong decay if they lie within 

1^ \ a few MeV of threshold. Furthermore, since the dominant component of 

C^ ■ these states is a superposition of two color-octet clusters which can be prO' 

O" 

f"^ I enhanced production of dibaryon resonances could be a signal of QGP for- 

mation. A total of eight, doubly-strange dibaryon states are considered for 



experimental search using the STAR detector (Solenoidal Tracker at RHIC) 
at the new Relativistic Heavy Ion Collider (RHIC). These states may decay 
to AA and/or p'E~ , depending on the resonance energy. STAR's large ac- 
ceptance, precision tracking and vertex reconstruction capabilities, and large 

^ , data volume capacity, make it an ideal instrument to use for such a search. 

^ I Detector performance and analysis sensitivity are studied as a function of 

resonance production rate and width for one particular dibaryon which can 
directly strong decay to pH~, but not AA. Results indicate that such reso- 
r> I nances may be discovered using STAR if the resonance production rates are 

C^ ■ comparable to coalescence model predictions for dibaryon bound states. 
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1 Introduction 

It is generally assumed that quantum chromodynamics (QCD) is capable of de- 
scribing the structure and spectroscopy of baryons and mesons. In addition, a 
number of models based on QCD predict other types of quark and gluonic systems 
such as pure gluonic states (glueballs), hybrids (qqg), four quark di-mesons (qqqq, 
qqqq) |jl|, penta-quark states (qqqqq), six quark dibaryons 0, ^, Q, strangelets, 
and for large systems several fm in size at sufficiently high temperature and/or 
net baryon density, the color-deconfined quark-gluon plasma (QGP). Extensive 
experimental searches have been made for all these systems. 

In particular, many experiments {e.g. |^, ^, ||, ^ 0) have searched for a stable 
(at least to strong decay), doubly-strange flavor- singlet dibaryon (Hq) of the six 
quark SU(3)-fiavor [SU(3)/] multiplets |^. This hypercharge (Y) = dibaryon 
was predicted by Jaffe to be stable against strong decay, but not to weak decay; 
hence the experiments generally searched for characteristic weak decay topologies. 
So far, no conclusive evidence for a bound Hq has been presented 0, §, |], ||, |^. 
However, the Hq mass may be larger than Maa = M\ + Ma, the AA strong decay 
threshold p, |10|, |TT[. In this case the Hq could be a strong interaction resonance. 



If the resonance lies only a few MeV above Maa, then it may be possible to see 
experimental evidence for its existence^ . In addition, it is important to point out 
that several other doubly-strange dibaryon states occur in the SU(3)/ multiplets 
which do not directly strong decay to AA and which may also appear as resonances 
in the AA and NE systems. 

Experimental searches for a stable Hq have attempted to create and observe 
the Ho through: (1) double-strangeness-exchange reactions, such as A{K~ ,K'^)X, 
A{K~ ,K^)X; (2) E^A capture p on nuclear targets with A>2; and (3) proton 
+ nucleus [^, 0] or nucleus + nucleus collisions where the Hq production mechanism 
could proceed via multi-hyperon production followed by coalescence |ll2|, ^ 0, |l^ 



Reactions such as d{K~ , i^'')Ho, for example, produce a relatively clean final state 
but suffer from poor momentum matching and subsequently, a reduced produc- 
tion rate. Heavy ion collisions, although producing a complex final state, gener- 
ate a large number of co-moving hyperons and non-strange baryons from which 
dibaryons could either coalesce or scatter through resonances. The invariant mass 
spectrum for AA pairs produced in 158 A GeV/c Pb + Pb collisions [data from 
WA97 ||TB|] does not reveal any resonance like structures. However, the data are 
too coarsely binned and lack sufficient statistics to be relevant for the type of ex- 
perimental search proposed here. On the other hand, a three standard deviation 
enhancement above background was observed in the AA invariant mass spectrum 

^ In the meson sector nature provides several examples of narrow states that lie a few MeV 
above their strong decay thresholds. The 0(1020), D*{2010)^, D*{2007)° and Dsi{2536)+ 
mesons are 32.1, 5.8, 7.1 and 35.0 MeV, respectively, above their strong decay thresholds and 
have total widths of 4.4, <0.13, <2.1 and <2.3 MeV, respectively. The latter three quantities 
represent upper limits at the 90% confidence level. 



from threshold to 30 MeV above threshold, in the AA pair production data from 
the ^^C{K-,K+) reaction [0. 

The Relativistic Heavy Ion Collider (RHIC) facility at the Brookhaven National 
Laboratory, together with four new detectors, is beginning a program of research 
whose primary goal is to produce a quark-gluon plasma and study its properties. 
It is noteworthy that a color-deconfined QGP could enhance significantly the pro- 
duction of Hq dibaryons compared to the production expected from a hot hadronic 



gas |jT3[ . The reason is that the largest component of the six quark dibaryon wave 
function consists of a color-singlet superposition of color-octet-octet components 

a ITS 




H) = J- \BB) + J- \8,®8,) . (1) 



Here, | 8c ® 8c) represents the color-octet-octet components, and | BB) denotes 
a summation of states composed of two physical baryons. Enhanced dibaryon 
production can occur by way of fusion of color-octet clusters within the deconfined 
plasma [0. If evidence for dibaryon resonances is observed in the data from 



RHIC collisions, then analyses of such data would provide an independent basis 
for studying the QGP. 

In this work we explore the feasibility of observing evidence for the existence 
of several doubly-strange dibaryon resonances through study of the AA and NS 
invariant mass spectra constructed from relativistic heavy ion collision data to be 
taken at RHIC.0 We focus on using the STAR |19[ detector [with the Silicon Vertex 
Tracker (SVT) and Silicon Strip Detector (SSD) [^ upgrades] since it is the most 
suited of the RHIC detectors for such a program. A preliminary study [^, pT[ | 



suggested that STAR could be used to search for Hq — > AA with resonance masses 
a few hundred MeV above Maa- The focus of the present work is on the NS decay 
channel with resonance masses a few MeV above the strong decay threshold. 

In Sec. H the doubly-strange dibaryons and their possible decay schemes are 
discussed. Using the P-matrix formalism the resonance widths are estimated 
in Sec. ^ In Sec. ^ the J'" = O"*", I(isospin) = 1, I3 = dibaryon member of the 
SU(3)/ 27-plet g [referred to as ^^(7^,1,13) with J^ = 0+, I = 1 and I3 = 0], 
which cannot directly strong decay to AA, is used to explore the statistical aspects 
and sensitivity of analysis of STAR data to the existence of such resonances, if 
they exist. A new, fast simulation detector response code (FSDR) is used for this 
study. Finally, a summary and conclusions are presented in Sec. 1^. 



^Formation of hypcr-nuclear baryon-baryon bound states or slightly unbound states is also 
possible {e.g. strange baryon analogs of the deuteron and the nucleon + nucleon ^Sq unbound 
state). Such states, if they exist, are expected to produce broad enhancements over many tens of 
MeV in the AA or iVS spectra, and not the narrow resonance structures predicted for six quark 
states which are the focus of the present work. 



2 Dibaryon States 

The direct product space of the J'^ = | baryon octet with itself can be written 
in terms of irreducible representations of SU(3)/ [H: 

8®8 = 1©8©8©10©T0©27. (2) 

The hypercharge [Y) ranges from +2 (NN, nucleon-nucleon) to —2 (SS) for these 
dibaryon states. The Y = 2 members include the deuteron of the 10-multiplet. 
The Y = 1 states contain NA and NS components. Sharp enhancements in the 
NA spectra at 2129 MeV fS^ [the NS threshold ^, |g] and 2139 MeV Q 
have been observed in single strangeness-exchange reactions such as d{K, 7r)X and 
d{TT,K)X. The Y= —1 and —2 states contain the experimentally difficult SA, SS 
and SS components. Finally, the Y = dibaryons contain AA, NS, AE, and SS 
components. Since A^s, As, Ss, and Hs are produced in relativistic heavy ion A + A 
central collisions, the RHIC experiments will provide an excellent opportunity to 
search for some of these Y = dibaryons. 

Table | lists the J'^ = 0+, 1+ y = dibaryon constituents of the various SU(3)/ 
multiplets in Eq. (^) from Refs. p|, || along with the dominant decay modes for 
a wide range of assumed dibaryon masses. The mass limits where strong and/or 
electromagnetic (EM) decay channels open are noted for each state. In general the 
best mass ranges to explore for resonances are those lying just above the strong 
and/or EM decay thresholds in Table |. Jaffe predicted a significant increase in 
the dibaryon masses with increasing dimensionality of the SU(3)/ representation. 
However, more recent calculations based on QCD sum rules I^S], ^ indicate that 
the Y = dibaryons in different multiplets should be similar in mass; the principal 
mass dependence of the 8®8 dibaryons is due to the explicit SU(3)/ symmetry 
breaking caused by the strange quark mass.0 All of the Y = dibaryon states 
in Table | will be considered in this work, since at present there is no reason to 
exclude them from experimental searches. 

For masses between M\\ and Mats the singlet Hq and H^^(0+, 0, 0) can strong 
decay to AA. Recently, using a P-matrix formalism, Bashinsky and Jaffe [Q consid- 
ered a hypothetical singlet Hq strong-decay resonance of mass several MeV above 
the A A (strong decay) threshold. For a variety of assumptions, they found that, 
owing to kinematic effects near threshold, the cross sections for S-wave AA elastic 
scattering should show structures whose scales are of the order of several MeV. 
It may be that the AA invariant mass spectrum generated from relativistic heavy 
ion collision data will show similar structure that can be taken as evidence of the 
existence of such a resonance. The analysis in @ also applies for the H^^(0+, 0, 0) 
case. 



^It is also worth noting that the Y — 2 deuteron member of the 10-flavor multiplet and the 
nucleon + nucleon ^Sq unbound state member of the 27- flavor multiplet differ in mass by only a 
few MeV. Calculations [pS based on a SU(3) invariant pseudoscalar and vector meson exchange 
model for the baryon-baryon system predict the 0+ 27-plet dibaryon mass to be less than the 
mass of the singlet Hq. 



The 8-, 10- and 10-plet J^ = 1+ dibaryons and the H27(0+, IJg) and H27(0+, 2,h 
dibaryons cannot strong decay to AA since antisymmetrization requires the fi- 
nal AA system to be spin-singlet-even or spin-triplet-odd and I(AA)=0. The 
H^'''(0+, 1, 0) and H^^(0"'", 2, 0) dibaryons can, however, strong decay to AA by way 
of the small 1=0 isospin admixture which contributes to the physical states due 
to Coulomb induced isospin mixing, but this transition rate should be much less 
than that for Hq —>■ AA. The four J'" = 1"*", l3=0 states can EM decay to AA7 
via El transitions leaving the AA in relative ^Pj states. The H^'^(0"'", 1, 0) and 
H^'^(0'*', 2, 0) states can also EM decay to AA7, for example via El transition to 
the ^Pi AA final state. All remaining states {i.e., those with I3 7^ 0) in this mass 
range (from Maa to M^-e) can decay only weakly. It is likely that the J'" = 1+ 
states are of larger mass than the 0"^ states and thus even more difficult to deal 
with experimentally. 

In the mass range from Mm- to M^s all states with I < 2 decay strongly to 
AA, NE or AE as indicated in Table | The H27(0+, 2, ±1) [H2^(0+, 2, 0)] states can 
strong decay to NE and AS [NE, AE and AA] via isospin admixtures or they can 
decay electromagnetically. The H^^(0"^, 2, ±2) remain blocked to both strong and 
EM decays up to M^s due to charge conservation. The H^^(0"'", 1, /s) members 
of the 27-multiplet cannot strong decay to AA because of isospin conservation; 
however, they can strong decay to NE and AE if the mass is greater than Mats 
and Mas, respectively. Above M^s all states in Table |I| can decay strongly. 

Table || lists the optimum decay channels and mass ranges for possible Y = 
dibaryon resonance searches. The last column indicates whether the channels 



include those appropriate for STAR [19|. Of the listed decay products, STAR 



can identify only protons (p). As and S~s. Eight of the Y = dibaryon states in 
Tables I and H (corresponding to each of the l3=0 states) are seen to be appropriate 
for a program of research with STAR. 

Evidence for the Hq and H^^(0"*", 0, 0) states might be seen in the AA invariant 
mass spectrum within several MeV of threshold. The four J"^ = 1"*", l3=0 states, if 
in the mass range from Maa to M„50, will decay electromagnetically to AA7. If the 
photon is not detected, the AA invariant mass spectrum would appear as a broad 
continuum. However, an enhancement still might be observed if the resonance 
masses are only a few MeV above Maa- These four resonances might also be 
observed in the pE~ spectrum if either lies just above the Mp^- = 2259.6 MeV 
strong breakup threshold. [] 

If the resonance mass exceeds Maa, the H^^(0+, 1, 0) and II^^(0+, 2, 0) states can 
decay either electromagnetically to AA7 or strongly via the small 1=0 admixture. 
The competition between these branches determines the experimental signature. 
If the strong interaction, isospin admixture decay dominates the EM decay, then 
the 1=1(2) state might be observed if its mass is between Maa and M„=o (M^oso). 

^Due to Coulomb interactions the nS'^ andpS~ mass thresholds are split by 5.1 MeV. Strong 
decay to NE therefore opens at M„20 =2254.5 MeV, whereas the observable decay channel for 
STAR, dibaryons pS", does not open until 2259.6 MeV. 



If the EM decay dominates, then both 1=1 and 2 states might only be observed 
as enhancements in the continuum just above Maa. Otherwise the H^'^(0^, 1,0) 
state might be seen in the pS~ invariant mass spectrum just above the Mp=- 
threshold. In this case H^^(0"'", 1,0) —>■ pE" might be observable with STAR. The 
H^'''(0"'", 2, 0) resonance could strong decay to TPTP and S+E~ final states, but the 
Ss, unfortunately, cannot be reconstructed by STAR. 

In this initial simulation study we prefer to consider the dibaryon resonance 
with smallest angular momentum which may strong decay to pS~ but not to 
AA, namely H^^(0"'", 1, 0). The reason for this choice is the expectation that the 
combinatoric background contribution to the pS~ invariant mass spectrum will be 
less than that for AA for the case of RHIC data for Au + Au central collisions (see 
Sec. Ij.) Simulation studies for the remaining dibaryon resonances in Table will 
be done in the future. 

It is worth noting that threshold-cusp effects |^, |2^, |2^ could produce en- 



hancements in the AA spectra at the A^H, SS, etc., thresholds, independent of any 
possible resonances. This is also true for the N'E spectra (1=1 component) at the 
AE and SS thresholds. Such enhancements are due to coupled-channels scattering 
effects near inelastic thresholds and must be treated carefully when searching for 
experimental evidence of dibaryon resonances. 

Resonances occuring between the nS° and pS~ thresholds will decay primarily 
to nS°, although decay to pH^ will occur if the resonance width is such that 
the distribution extends above Mps-- For resonance masses near or a few MeV 
above Mp=- , the branching ratio between the nS° and pH~ channels depends upon 
the dynamics of the raS° — pS~— dibaryon coupled channels system, as well as on 
the kinematic density of final states factor. For resonance energies much above 
MpH- the branching ratio is determined mainly by isospin invariance which results 
in a 50% decay fraction for both channels. Thus the study in Sec. ^ will be 
in terms of the dibaryon resonance width and production rate for the observable 
dibaryon— >■ pS~ decay channel only. Inclusion of the model-dependences associated 
with the nS° — pS~— dibaryon system lies well beyond the scope and purpose of 
this work. 

Finally, we note that the decay products of resonance states formed deep within 
the interior of the collision region will undergo strong, final state interactions (FSI) 
so that observation of these resonances becomes more difficult. It is likely that only 
resonances occuring in the periphery of the QGP or in the mixed QGP-hadronic 
phase prior to freeze-out will lead to experimental signatures in invariant mass 
spectra. Detailed calculations of FSI effects on this and other potential QGP 
signals remain to be done. 

3 Dibaryon Resonance Widths 

Bashinsky and Jaffe used the P-matrix scattering formalism to compute the 
AA elastic scattering S-wave total cross section in the presence of a flavor-singlet, 
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Hq dibaryon resonance located a few MeV above the AA threshold. Their model 
is used in this section for the similar case of a possible H^^(0+, 1,0) (for brevity 
referred to as H^^ in the remainder of this paper) dibaryon resonance lying just 
above the iVS threshold. 

The total cross section for iVS elastic scattering is given by 

o-TOT = ^-Ini(/ii) , (3) 

where /n is approximated by the forward, S-wave elastic scattering amplitude, 
which in the P-matrix formalism is written as 

"•I 

The reduced K-matrix from Ref. P] as a function of total energy (e) in the NE 
channel is 

Kr\^)^-K — ^1^^ . (5) 

e — ^thi — ^1 
where, 

_, fli , 1 - Pibo + 2meiaibo , ci 

1 - pibo (1 - piboy 1 - pibo 

In Eqs. (§)-(§), ki is the relative center-of- momentum system wave number, m 
is the reduced mass of the two-baryon (NE) system, ethi is the NS threshold 
energy, ei = e.^ — ethi, and e^ is the assumed resonance mass of the H^^ dibaryon. 
Parameter b^ is a resonance size parameter with typical value (150 MeV)~^, pi is 
the resonance width parameter which is a function of the resonance state - decay 
channel coupling, and ai is the NE scattering length which was assumed to be 
(200 MeV)-^ 

Elastic total cross section results obtained using Eqs. (H)-(|l) are shown in Fig. ^ 
by the solid lines for assumed values of ei = 2, 3 and 5 MeV and for the typical 
resonance width parameter pi=50 MeV. The unitarity limit, Air/kf, is indicated by 
the dashed lines. These results indicate that the resonance distribution displays a 
characteristic width roughly proportional to the difference between the resonance 
and threshold energies. 

Resonances decaying to pS~ may also decay to nS° which is 5.1 MeV lower in 
energy. These results suggest that resonances near Mp^- would be about 5 MeV 
wide. It is reasonable to assume, however, that the pS~ spectra for such cases 
would display a narrower range of enhancement because of the kinematic limit 
which cuts off the lower energy portion of the pS~ invariant mass distribution. 



4 Search for Dibaryon Resonances using STAR 

An environment conducive to the formation of dibaryon resonances may be the 
hot, dense nuclear matter that will be created through central Au + Au collisions 
at RHIC The experimental objective, then, is to infer evidence for the existence of 
such resonances through examination of AA or pE~ invariant mass spectra. Since 
the resonance cross sections are expected to be small, a large acceptance detector 
with the ability to do precision tracking to within a few cm of the production 
vertex is needed. In addition, the detector must handle high data rates and have 
event processing capabilities for very large data volumes. The STAR detector 



| T9[] at RHIC with its central vertex tracking system (SVT-SSD) is well-suited for 
such experiments. For the simulations to be discussed here, the parameters in 
Table fT^ were used to approximate the kinematic acceptance and reconstruction 
performance of STAR (including the SVT-SSD). The momentum resolution and 
track reconstruction efficiencies were taken from previous simulation studies [ ppj . 
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Coalescence model calculations |13, 15| for Au-(-Au central collisions at AGS 



energies (~11 A-GeV fixed target) predict a wide range of production rates for the 
fiavor-singlet Hq from of order 10~^ to 0.5 per event, depending on the assumed 
temperature and strangeness content of the hadron gas and on the assumed AA 
interaction model. Color- deconfinement in the QGP could enhance the production 
rate via the color-octet-octet component of the dibaryon wave function [Eq. (Q)]. 
At RHIC energies increased strangeness production and possible QGP formation 
could further enhance the dibaryon production rate above this range. Source 
expansion and collective fiow in RHIC collisions could also affect coalescence model 



predictions []14[. Comprehensive calculations including all these effects have not 
been done. Nevertheless, the results of such calculations suggest the necessary 
level of sensitivity for the experimental search to be meaningful. 

In the following, an analytical model is presented which relates the statistical 
significance of the dibaryon resonance signal to the number of Au+Au central 
RHIC collision events analyzed and the assumed resonance production rate and 
width. This is followed by a description of a fast numerical simulation model and 
presentation of statistical and sensitivity results for experimental detection of a 
hypothetical, narrow dibaryon resonance in the pH" channel just above threshold. 

4.1 Analytical Model 

The integrated signal for the H^^ — > pS^ invariant mass distribution can be ap- 
proximated by 

S = Y^ n'^epesfn-cut = rinepe-EfH-cutNev = SN^v , (7) 

1=1 



where n^ is the assumed number of dibaryon resonances in the i^^ Au+Au central 
colhsion event that decay to pE~ with all daughter products going fully into the 
kinematic acceptance range specified in Table ^, and uh is the average number 
of H^^ — >■ pS~ decays per event for the data sample consisting of Nev events. Also 
in Eq. (|^ Cp and 6= are the proton and S~ reconstruction efficiencies, respectively 
(number correctly reconstructed and accepted divided by number in the kinematic 
acceptance or detector fiducial volume), fn-cut is the fraction of reconstructed 
jj27 _^ p'^- resonance decays that survive any remaining analysis cuts {e.g. cuts to 
remove fiuctuations from the background subtracted spectra and the upper/lower 
pS~ invariant mass limits), and S is the average signal per event |2^. The total 
background underneath the signal peak is approximated by 

B = Y.iK^^ - ^HepeE)F{pi) = BN,, , (8) 

where iV^ and A^^ represent the number of reconstructed and accepted proton and 
S~ candidates for the i^^ event, F{pi) represents the relative fraction of random 
pS~ pairs whose invariant mass occurs within the domain of the Yi^^ mass peak, and 
B is defined as the average background per event. The fraction, -F(pi), depends 
on the proton and S~ momentum space distributions for the collision events (after 
cuts) and on the pS~ invariant mass range which in turn depends on the width of 
the resonance peak. 

The numbers of candidate protons and S" include several contributions which 
are represented by the following: 

N; = {n'j^ + u'h^p + nif^^p)ep + NFaise,p (9) 

m = {n'H + n^H,E+^lh,E)e^ + NFalse,E, (10) 

where n^ (^ks) is the number of protons (S~s) for the i*'* event produced by 
H^^ resonance decays for which the proton (all decay products of the H~) entered 
the detector acceptance but all decay products of the E~ (the proton) did not. In 
Eqs. (ID and (|1^) nlf^^ and nj^ = represent the remaining number of protons and 
S~s produced in the collision, for example by hadronization and/or rescattering 
processes, which enter the detector acceptance. The quantities NFaise,p and Npaise,- 
represent the number of incorrectly identified and accepted protons and S~s per 
event. The numbers of proton and H~ contaminants depend on various analysis 
cuts, overall event multiplicity, etc. In general, tighter cuts in the reconstruction 
and analysis result in lower efficiency and reduced contamination, while relaxed 
cuts have the opposite effects. The efficiencies and contaminants per event in this 
section are averages of inclusive quantities. 

The total number of counts, or yield {Y), for the invariant mass spectra in the 
domain of the H^^ resonance peak for Nev Au+Au collision events is 

Y = ST + B = S + S' + B = {S + S' + B)Ney, (11) 



where 

'S't = X! ^^H^p^s — ^H^-pe-Nev = S + S' (12) 

j=l 

which defines S'. The total error in the signal is assumed to be a sum of statistical 

and systematic errors given by 

AS = ASstat + ASsyst , (13) 

where the statistical error in 5* is determined by the independent statistical errors 
in the total yield (AY), background (AB), and S' {AS'), where 



ASstat = ^{AY? + {AS'f + (Ai?)2 = ^5 + 2(5 + S') . (14) 

The systematic error in S is assumed to be proportional to A^e^,, where ASsyst = 

^>^syst-^^ ev 

The statistical significance of the observed resonance peak can be expressed as 
the ratio of signal to signal error, which in the present model is given by 

o ST\g^^ , , 

(15) 



^S Js + 2{B + S')^/NZ + ASsystNev 



If fu-cut is constant [^ and systematic errors are negligible, S/ AS increases with 
\/Nev as expected, whereas for finite ASsyst, S/AS is limited to S/ ASsyst for many 
events. 

Solving Eq. (|l^) for N^v yields (assuming constant fu-cut HH) 

_ {S/ASnS + 2{B + S')] 
^^ [S- {S/AS) ASsyst? ^ ^ 

which provides an estimate of the number of collision events which must be ob- 
tained and analyzed in order to achieve a specified value of S/AS. Systematic 
errors limit the detection capability such that 

S > {S/AS)ASsyst , (17) 

which results in a minimum, detectable dibaryon resonance production rate given 
by, 

n„ > ^Il^^j^i^, (18) 

assuming the required number of events from Eq. (|I6[) can be achieved. 

Optimization of a particular experiment and data analysis program {e.g. num- 
ber of events and cut parameters) can be guided by the results presented in this 
section provided a reasonable estimate of systematic errors can be made. The 
latter must wait until experience is gained with analysis of actual STAR data and 
therefore lies well beyond the scope of the present study. The remaining discussion 
and results focus on the statistical requirements for detecting H^^ -^ pS~. In the 
following, Eqs. (piSJ) and ([T^) were used to estimate the dependence of {S/AS) on 
Nev and uh- 

10 



4.2 Numerical Simulations 

A fast simulation detector response code (FSDR) was developed to provide rapid 
evaluation of the capabilities of relativistic heavy ion experiments with respect to 
proposed new measurements and research programs. In simplest terms, FSDR 
projects a given input list of particles from an event generator directly to the final, 
reconstructed particle list. Specific detector acceptance and track reconstruction 
performance parameters are supplied from realistic simulations and analyses. The 
values assumed for STAR are hsted in Table fT| |^, ^ Q. FSDR propagates 



(freely) and decays unstable particles such as A and H~ according to known branch- 
ing ratios and lifetimes and applies acceptance cuts, track finding efficiencies, par- 
ticle identification efficiencies, and momentum resolution smearing to the charged 
primary particles and to the charged daughter particles from the decays. 

Following this, FSDR reconstruction of A and H~ is very similar to that used 
in the actual STAR reconstruction analysis. The charged particles selected for this 
analysis were assumed to be reconstructed in both the SVT-SSD and TPC (Time 
Projection Chamber). The SVT-SSD provides excellent track position resolution 
(few tens of microns) [pO| , p8| , |30| in the decay region within several cm of the 



primary vertex. Reconstructed trajectories (assumed to be helices) were projected 
to the primary vertex where an impact parameter cut (1 mm) differentiated pri- 
mary particles {i.e., those assumed to emerge from the primary collision vertex) 
from secondary particles {i.e., those assumed to come from decay vertices). Can- 
didate A decays were found by a distance of closest approach (DCA) cut (2 mm) 
among the projected trajectories of the secondary protons and 7r~s. The result- 
ing A reconstruction in the pTT~ invariant mass spectrum is shown in Fig. ^ for 
600 central Au -|- Au hijet [^ collision events which included an average of 6.6 
jj27 _^ p^- __^ p{A7T~) embedded decays per event (see following discussion). The 
cross-hatched region indicates the number of reconstructed p — 7c~ pairs actually 
produced by A decays in the simulation. The reconstruction efficiency for pri- 
mary As {i.e. number correctly reconstructed divided by number in acceptance) 
obtained here was about 5-6%,^ in reasonable agreement with the 8% value in 



Ref. [0 and the range of 4 — 15% in Ref. ||20[. These comparisons verify that 
the FSDR A simulation results presented here are consistent with those obtained 
from more detailed studies. The ±3cr width of the reconstructed A peak in Fig. ^ 
is about 8 MeV which quantitatively agrees with that found in Ref. ||30|. Further 
cuts could reduce the background for primary As as in Ref. [Q, but this was not 
done here in order to retain the secondary As from S~ decays. 

A DCA cut (4 mm) among all pairs of selected A candidates [i.e., those with 
impact parameter less than 2 cm from the primary vertex and which have a re- 
constructed mass within ±3cr (±4 MeV) of the peak in Fig. Q with all remaining 

®The reconstruction efficiency for As from E!~ decays is somewhat higher than that for primary 
As since the decay vertices are distributed farther from the primary vertex where the proton and 
TT^ daughters are not as hkely to be removed by the 1 mm impact parameter cut as occurs for 
the daughter tracks from primary A decays. 
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secondary 7r~s yielded the S~ spectrum shown in Fig. ^ using the same events as in 
Fig. 1^. The cross-hatched region indicates A and vr" reconstructed pairs that orig- 
inate from S~ decays. The FSDR S~ reconstruction efficiency with the preceding 
cuts was about 0.06 which is larger than that in Ref. [Q. The higher reconstruc- 
tion efficiency obtained here is offset by the smaller peak signal-to-background 
ratio here (~ ^) compared to the larger value (~ 1) in Ref. |30]. The ±3o" width 



of the reconstructed H peak in Fig. ^ is about 10 MeV, in good agreement with 



Ref. IQ. We refer to the present set of cuts as "relaxed" while those in Ref. [0 
are "tight." A tight selection of primary protons (impact parameter from primary 
vertex < lOO/im) and S~ candidates (within ±6 MeV of the peak in Fig. |^) with 
projected DC As < 6 mm, provided the set of pS~ pairs used in the invariant mass 
plot in Fig. |. 

Variable numbers of H^^ resonances were randomly embedded in the hijet [^ 



event generator output for 200 A-GeV Au-|-Au central collisions according to the 
following distribution, 



In Eq. (|T9]) m^ = yjM{H'^'^y + p^^ is the transverse mass, T = 238 MeV, p = 
3.0 GeV~^, y4 is a normalization constant, and at midrapidity {y=0) the rriT distri- 
bution corresponds to an effective temperature of 139 MeV. For this distribution 
function approximately 39% of the total H^^s have all four decay particles fully 
contained within the acceptance (Table |IT|). In the following discussion, produc- 
tion rates refer to the number of H^^ -^ pE^ decays per event for which all four 
decay particles are within the STAR acceptance. 

In FSDR the H^'^ masses were randomly distributed according to the P-matrix 
resonance mass distribution (Sec. |^). The latter was obtained from the NE total 
elastic cross section in Eq. (^ by eliminating the non-resonant amplitude (ai was 
set to 0) and by removing the incident flux factor in the definition of the total 
cross section {i.e. by multiplying by ki/m). The resonant pS" invariant mass 
probability distribution was therefore assumed to be 

P(Mh27) = PoA:i^TOT(e, a, = 0)/m , (20) 

where Mh27 = e is the randomly sampled H^^ mass and Pq is a normalization 
constant. 

In this work the H^^ resonance energy, ei, was assumed to be 2 MeV above the 
pS~ threshold and a range of width parameters, pi, was assumed where pi = 10, 
23, 37 and 50 MeV, corresponding to mass distributions with full width at half 
maximum (FWHM) values of 0.94, 2.4, 4.1 and 5.6 MeV, respectively. Average 
H^^ production rates (nn) of 3.1 (8 total) and 4.7 (12 total) per central Au+Au 
HIJET event were assumed. For each of the resulting eight cases, which correspond 
to different values for pi and hh, 1940 events were analyzed. 

The preceding H^^ production rates and number of events were selected in or- 
der to generate statistically significant resonance signals, with modest computing 
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requirements, such that straightforward background subtraction and signal deter- 
mination techniques would suffice. The analytical model in the preceding subsec- 
tion, when supplemented with the present FSDR results, can be used to obtain 
estimates of STAR's detection sensitivity and data volume requirements for much 
smaller dibaryon resonance production rates. 

The pE~ invariant mass spectrum is shown in Fig. ^ for 1940 Au+Au central 
collision events for the case in which pi = 23 MeV (FWHM = 2.4 MeV) and 
uh = 3.1. The H^^ peak is quite apparent for this case when the resonance occurs 
near threshold and the background is rapidly decreasing. The contributions of 
the large number of false S~s (see Fig. H) are dispersed throughout the spectrum. 
Notice that momentum resolution effects do not degrade the H.^^ peak. Finite 
momentum resolution |2^ for the final, reconstructed daughter protons and pions 
in the decay chain, dibaryon -^ pE~ — >• p{An~) -^ p{pn~)n~, results in about 
10 MeV/c uncertainty in the pS~ relative momentum. For this case, where FWHM 
is 2.4 MeV, such effects would only broaden the peak in the reconstructed p'E~ 
invariant mass spectrum to about 2.6 - 2.9 MeV. 

The uncorrelated pS~ background distribution near threshold is proportional 
to y/e — ethi- Therefore, we assumed the threshold constrained model for the pS~ 
background distribution given by 

Bmod{x) = aiy/x + a2x'^ + a^x^ + a^x^ + a^x^, (21) 

where x = t — ethi and Oi, 02, ... 05 are parameters to be determined by chi-square 
fits. Excellent fits were obtained to the nonresonant backgrounds for all cases as 
exemplified by the fit shown in Fig. | by the solid line. The cross hatched portion 
of the spectrum in Fig. § indicates the reconstructed pS~ pairs from actual H^^ 
decays where the remaining background is seen to be in quantitative agreement 
with the fit. 

The background subtracted peak is shown in the inset panel in Fig. ^ For this 
and all other cases the mass range of the resonance peak was taken from threshold 
to 7 MeV above threshold. In addition, a lower cut-off of 20 counts per 1 MeV 
bin (specific for the number of events used in this analysis) was applied to remove 
the residual background fluctuations which are apparent in the inset panel. The 
remaining counts in the 7 MeV invariant mass range, which were above the 20 
counts per bin cut-off, constitute the measured signal, Sfsdr- For this case the 
H^^ reconstruction efficiency is 0.047. 

The results for each case are summarized in Table |I^ and in Fig. |^ which give 



the statistical significance of the signal, {S/AS)fsdr- The quantity {S/AS)fsdr 
was calculated using 

f S \ Sfsdr . 



\ASJ FSDR ^J2Yfsdr - Sfsdr ' 

where Yfsdr is the total yield in the 7 MeV wide domain of the resonance peak. 
Also listed in Table |^ are estimates of the number of Au+Au events which are 
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needed in order to achieve {S/AS)fsdr = 3, assuming statistical errors only, 
constant fn-cut, and using the N^v dependence from Eq. (|l^) to scale the FSDR 
results. 

The decrease in {S/AS)fsdr with increased resonance width and fixed pro- 
duction rate is due to increased signal losses {i.e. reduced fn-cut) where, (1) more 
counts are lost in the tails of the resonance distribution above the upper mass limit, 
and (2) more counts are lost due to the 20 counts per bin cut-off. If the resonance 
mass domain was also increased, then {S/AS)fsdr would further decrease due to 
the larger background contribution to AS*. From Table ^ and Fig. ^ it is seen that 
the required number of events to analyze in order to achieve a certain statistically 
significant signal [e.g. {S/AS)fsdr = 3] increases for broader resonances and/or 
reduced production rates, as expected. 

For analysis of real STAR data, it is likely that a more accurate background 
subtraction method will be needed than that used here, since optimistic production 
rates were assumed in this numerical study. A possible method is to form random 
pS~ invariant mass histograms using mixed event pairs, normalized to the actual 
number of pS~ pairs in the data. This removes dynamical correlations in the pS~ 
background spectra and is analogous to similar methods used successfully in pion 
interferometry analyses [^ 



For smaller rates of production {e.g. uh ~ 1) Eq. ([T6|) simplifies to, 

nr ^ 2{S/ASfnth,pep{nth,EeE + NFaise,E)F{pi) 
[nHepe-tH-cut] 

where contributions from ASsyst and S' were neglected. The required number of 
events is proportional to n']^ as expected [^. For n/f(total)~l and using the 



results in Table |V|, it is estimated that of order 10^ — 10^ Au+Au events must 
be analyzed to achieve a minimum dibaryon signal with S/ AS = 3. At STAR, 
nominal run-time operations should obtain 10^ central Au-|-Au events in about 
two weeks. 

It is likely that the more crucial limiting factor in STAR's ability to detect 
dibaryon resonances, if they exist in the mass range considered here, will be due to 
systematic errors. Although the sources and magnitudes of these types of errors are 
not known at this time, a simple estimate can be made assuming the false counts 
due to systematic errors scale with overall yield. From Eq. (|T^ the minimum 
detectable signal is determined by the inequality in Eq. (|1^ provided N^^ can be 
arbitrarily large. If we assume ASgyst = fB, where / is a constant, and use Eq. ([^ 



with constant JH-cut I29| , the detectable range of dibaryon resonance production 
is given by 

UH > {S/AS)fnH{B/S) (24) 

where hn and S refer to a specific simulation result. The n/^ (total) =8 FSDR 
results indicate that {B /S) ~ 3. Requiring S/ AS = 3 and assuming, for example, 
that the systematic errors are 1% of the background, the minimum, detectable 
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dibaryon production rate is n// (total) ~ 0.7. Since STAR is a dedicated, long-term 
facility, discovery of dibaryon resonances of the type discussed here is well within 
reach of the experiment, provided the systematic errors can be kept small (to a 
few percent) and the dibaryon production rate is of order 1 per event or greater. 

5 Summary and Conclusions 

The possibility that Y = dibaryon states may occur as narrow resonances, if 
located in energy just above the strong interaction breakup threshold, offers inter- 
esting new discovery opportunities for the relativistic heavy ion physics program 
at RHIC, especially for the STAR experiment. In the baryon octet octet direct 
product space for F = dibaryons we have hsted eight possible states (correspond- 
ing to 13=0) which could, in principle, be discovered by STAR with its SVT-SSD 
central vertex tracking system. If produced, these l3=0 states decay to AA and/or 
pS~ channels, depending on the resonance energy. In addition, the formation of 
these resonances may be significantly enhanced in a color-deconfined medium, such 
as the QGP, by way of the dominant color-octet-octet component of the dibaryon 
wave function. A significant change in the production rate of these states, in 
conjunction with other QGP signals, would provide strong corroborative evidence 
for QGP formation. High quality statistics for the pS~ invariant mass spectrum 
can be expected at STAR with the analysis of data from of order several hundred 
thousand central Au-|-Au collision events — several days of STAR data acquisition. 

One specific Y = dibaryon state was used to illustrate and estimate the 
detection and sensitivity range of STAR. This state, the J'^ = 0"*", 1=1, l3=0 
member of the 27-plet, was assumed to lie within a few MeV above the pS~ 
strong breakup threshold. The resonance distribution was estimated using the P- 
matrix formalism 0, and the detector response and event reconstruction/analysis 
effects were estimated using a fast simulation, detector response code. This new 
simulation tool, FSDR, enables rapid evaluation of relativistic heavy ion detec- 
tor/reconstruction/analysis capabilities for proposed, new physics programs at 
RHIC. 

The ability of the STAR experiment to detect a possible H^^ -^ pS^ dibaryon 
resonance decay was evaluated in terms of the statistical significance of the dibaryon 
signal as a function of the H^^ production rate and width (FWHM). The numerical 
simulation results, together with an analytical model, were used to estimate the 
data volume requirements for much smaller dibaryon resonance production rates. 
Improved background subtraction methods were suggested for the case of reduced 
production rates; limitations due to systematic errors were also estimated. The 
large data volume required (of order 10® central events) for rare particle searches 
of the type discussed here is compatible with the nominal run plan for STAR. The 
discovery of possible dibaryon resonances with widths of order a few MeV which lie 
just above their strong decay threshold is feasible at STAR, provided systematic 
errors are small (few percent) and the dibaryon resonance production rate is of 
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order 1 or greater per Au+Au central collision. 
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Table I: J'^ = 0+, 1^ Y = dibaryon constituents of the 808 baryon octet - 
baryon octet direct product space for each irreducible SU(3) flavor representation. 
The dominant decay modes for each resonance mass range are shown. NN ab- 
breviates "mass of nucleon + nucleon," and similarly for the other baryon-baryon 
pairs. Weak decays which reduce the number of strange quarks by 1 or 2 are de- 
noted by "AS=1" or "AS=2," respectively. "S" denotes strong interaction decays. 
"EM" denotes electromagnetic decays. "S-Iso" indicates strong interaction decays 
via the small isospin admixtures in the physical states. "S"' indicates that strong 
interaction decays from the 1=1 states are possible to both NS and AS. 



SU(3)^ 
Irrep. 


(^^I,I3) 






Resonance Mass 


Range 






NN-NA 


NA-NS 


NS-AA 


AA-iVS 


iVS-AE 


AS-SS 


>EE 


1 


0+0 


AS=2 


AS=1 


AS=1 


S 


S 


s 


S 


8 


1+0 


AS=2 


AS=1 


AS=1 


EM 


s 


s 


S 


8 


1+ 1 1 


AS=2 


AS=1 


AS=1 


AS=1 


s 


s 


s 


8 


1+ 1 


AS=2 


AS=1 


AS=1 


EM 


s 


s 


s 


8 


1+ 1 -1 


AS=2 


AS=1 


AS=1 


AS=1 


s 


s 


s 


10 


1+ 1 1 


AS=2 


AS=1 


AS=1 


AS=1 


s 


S' 


s 


10 


1+ 1 


AS=2 


AS=1 


AS=1 


EM 


s 


S' 


s 


10 


1+ 1 -1 


AS=2 


AS=1 


AS=1 


AS=1 


s 


S' 


s 


10 


1+ 1 1 


AS=2 


AS=1 


AS=1 


AS=1 


s 


S' 


s 


10 


1+ 1 


AS=2 


AS=1 


AS=1 


EM 


s 


S' 


s 


10 


1+ 1 -1 
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AS=1 


AS=1 


AS=1 


s 


S' 


s 


27 


0+0 
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AS=1 


AS=1 


S 


s 


s 


s 


27 


0+ 1 1 


AS=2 


AS=1 
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AS=1 


s 


S' 


s 


27 


0+ 1 


AS=2 


AS=1 


AS=1 


EM, 

S — Iso 


s 


S' 


s 


27 


0+ 1 -1 


AS=2 


AS=1 


AS=1 


AS=1 


s 


S' 


s 


27 


0+2 2 


AS=2 


AS=1 


AS=1 


AS=1 


AS=1 


AS=1 


s 


27 


0+ 2 1 
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AS=1 


AS=1 


AS=1 


EM, 

S-Iso 


EM, 
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S — Iso 
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EM, 
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AS=1 
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Table II: Optimum decay channels and resonance mass ranges for J'^ = 0^,1^ Y 
dibaryon resonance searches. 



SU(3); 






Accessible 


Irrep. 


iJ\l,h) 


Decay Channels and Resonance Mass Ranges (MeV) 


to STAR 


1 


0+0 


AA ~ 2231 


Yes 


8 


1+0 


AA ~ 2231; nE° ~ 2254; pE' ~ 2260 


Yes 


8 


1+ 1 1 


pS° ~ 2253 




8 


1+ 1 


AA ~ 2231; nE^ ~ 2254; pE~ ~ 2260 


Yes 


8 


1+ 1 -1 


nE' ~ 2261 




10 


1+ 1 1 


pS° ~ 2253 




10 


1+ 1 


AA ~ 2231; nE° ~ 2254; pE' ~ 2260 


Yes 


10 


1+ 1 -1 


nE- ~ 2261 




10 


1+ 1 1 


pS° ~ 2253 




10 


1+ 1 


AA ~ 2231; nS° ~ 2254; pE' ~ 2260 


Yes 


10 


1+ 1-1 


nE- ~ 2261 




27 


0+0 


AA ~ 2231 


Yes 


27 


0+ 1 1 


pH° ~ 2253 




27 


0+ 1 


AA from 2231 to 2254;gnS° ~ 2254; pE' ~ 2260 


Yes 


27 


0+ 1 -1 


nE' ~ 2261 




27 


0+2 2 


S+S+ ~ 2379 




27 


0+ 2 1 


pS° from 2253 to 2382;g S+S^ ~ 2382 




27 


0+2 


AA from 2231 to 2385;" SS ~ 2385 


YesQ 


27 


0+ 2 -1 


nE' from 2261 to 2390;g E^E' ~ 2390 




27 


0+ 2 -2 


S-S- ~ 2395 





° If strong decay via isospin admixture dominates EM decay and resonance remains narrow; 
otherwise AA for resonance mass ~ 2231 MeV only. 

'' If strong decay via isospin admixture dominates EM decay and resonance remains narrow; 
otherwise p2° for resonance mass ~ 2253 MeV only. 

■^AA decay channel only. 

'^ If strong decay via isospin admixture dominates EM decay and resonance remains narrow; 
otherwise nE' for resonance mass ~ 2261 MeV only. 



21 



Table III: Kinematic acceptance and reconstruction parameters assumed in the 
analytical and numerical simulations. 



Parameter 


Value 


Transverse momentum acceptance 


Pt > 0.1 GeV/c 


Pseudorapidity acceptance 


\r] <1 


Azimutlial acceptance 


(j) = 0^2n 


Momentum resolution, Ap/p 


2% 


TPC track reconstruction efficiency^ 


90% 


SVT track reconstruction efficiencyf] 


80% 


Particle Identification (tt, protons) 


ioo%0 



"Time Projection Chamber; main tracking detector for STAR. 

''Silicon Vertex Tracker and Silicon Strip Detector tracking system for STAR. 

■^Actual value not yet available. 
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Table IV: FSDR results for the statistical significance {S/AS)fsdr of the H^'^ -^ 
pE~ dibaryon resonance decay signal for total production rates per Au+Au central 
event of 8 and 12 and for four assumed values of resonance width. 



Uh 


nHf\ 


Pi 


FWHM 


{S//\S)fsdi^ 


Nev for 




(Total) 




(MeV) 


(MeV) 






{S/AS)fsdr -- 


= 3E 


8 


3.1 


10 


0.94 


7.9 




280 




8 


3.1 


23 


2.4 


6.4 




426 




8 


3.1 


37 


4.1 


5.6 




557 




8 


3.1 


50 


5.6 


4.1 




1039 




12 


4.7 


10 


0.94 


10.2 




168 




12 


4.7 


23 


2.4 


9.1 




211 




12 


4.7 


37 


4.1 


8.6 




236 




12 


4.7 


50 


5.6 


7.5 




310 





'^Number of H^^ -^ pS^ decays per event in which aU decay products enter the STAR detector 
acceptance. This is 39% of the total for the assumed H^^ momentum distribution model. 
''For 1940 Au+Au central collision events. 
■^Assuming {S/AS)fsdr. is proportional to \/Nev as discussed in the text. 
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Figure 1: Total elastic cross sections for A^+S scattering assuming the H^^ dibaryon 
resonance is 2, 3 and 5 MeV above threshold are shown in the left, center and right 
panels, respectively. The dashed lines indicate the unitarity limit, An /k\. The N'E. 
threshold energy is denoted by e^p (axis label). 
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Figure 2: FSDR p7r~ invariant mass spectrum for A -^ pn" reconstruction for 600 
central Au+Au hijet events with embedded H^^ resonances at a collision energy 
of 200 A-GeV as described in the text. The cross-hatched portion indicates the 
number of correctly reconstructed A -^ p7r~ decays. Bin size is 2 MeV. 
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Figure 3: FSDR Att^ invariant mass spectrum for S^ ^ Att^ reconstruction for 
600 central Au+Au hijet events with embedded H^^ resonances at a collision 
energy of 200 A- GeV, using the relaxed S^ decay vertex reconstruction cuts de- 
scribed in the text. The cross-hatched portion indicates the number of correctly 
reconstructed S~ -^ Avr" decays. Bin size is 2 MeV. Note the suppressed zero for 
the vertical scale. 
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Figure 4: FSDR pS~ invariant mass spectrum for H^^ -^ pS^ reconstruction for 
1940 central Au+Au hijet events at a collision energy of 200 A-GeV, assuming 
an average of 3.1 H^'^ -^ pS~ resonance decays per event into the acceptance with 
resonance energy 2 MeV above threshold and P-matrix width parameter (pi) of 
23 MeV (FWHM = 2.4 MeV). The cross-hatched portion indicates the number of 
correctly reconstructed H^^ — > -pEr decays. The background subtracted peak is 
shown in the inset panel. Bin size is 1 MeV. Solid curve indicates the five-term 
threshold constrained background model fit. 
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Figure 5: Statistical significance of the H^^ — ^ pH" dibaryon resonance decay signal 
as a function of resonance width (FWHM) for resonance production rates of 3.1 
(triangles) and 4.7 (dots) per event (all decay products enter detector acceptance). 
For each case 1940 Au+Au central collisions were analyzed. Solid lines connect 
the values for fixed production rates. 
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